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In the standard scenario of planet formation, planets are formed from a protoplanetary 
disk that consists of gas and dust. The building blocks of solid planets are called as planetes- 
imals that are formed by coagulation of dust. We review the basic dynamics and accretion 
of planetesimals by showing TV-body simulations. The orbits of planetesimals evolve through 
two-body gravitational relaxation: viscous stirring increases the random velocity and dynam- 
ical friction realizes the equiparation of the random energy. In the early stage of planetesimal 
accretion the growth mode of planetesimals is runaway growth where larger planetesimals 
grow faster than smaller ones. When a protoplanet (runaway-growing planetesimal) exceeds 
a critical mass the growth mode shifts to oligarchic growth where similar-sized protoplanets 
grow keeping a certain orbital separation. The final stage of terrestrial planet formation is 
collision among protoplanets known as giant impacts. We also summarize the dynamical 
effects of disk gas on planets and the core accretion model for formation of gas giants and 
discuss the diversity of planetary systems. 
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§1. Introduction 

The solar system consists of planets, their satellites, and a huge number of 
minor bodies (asteroids, trans-Neptunian objects, and comets). The planets can be 
classified into three groups: terrestrial planets (Mercury, Venus, Earth and Mars), 
gas giants (Jupiter and Saturn), and ice giants (Uranus and Neptune). These groups 
differ from one another by compositions, planetary masses, and orbital radii. The 
terrestrial planets are light (masses M < M^, where is the Earth mass, 6 x 
10^'' g), rocky ones with relatively small orbital radii (a ^ 1 AU, where AU is 
the astronomical unit that is the average distance between the sun and the Earth, 
1.5 X 10^'^ cm), gas giants are heavy (M > 10^ M^) planets with main components of 
H/He gas and a ~ 5-10 AU, and the ice giants are moderately massive (M ~ lOM^ 
with main components of H2O/CH4/NH3 ice and are in distant regions (a ~ 20-30 
AU). These planetary orbits are nearly circular and coplanar, which suggests that 
the solar system was formed from a protoplanetary disk around the proto-sun. 

The standard scenario for solar system formation was established in 1960's to 
igSO's.S'l^ The basic ideas of the standard scenario are: 

1) disk hypothesis: A planetary system is formed from a protoplanetary disk with 
mass much less than its host star's mass, 

2) planetesimal hypothesis: Building blocks of solid planets are km-sized rocky/icy 
"planetesimals" that are formed by coagulation of dust grains in the disk, 

3) core accretion model: Gas components of gas giants are added after rocky/icy 
planets (cores) accrete from the planetesimals. 

Though the standard scenario still has serious difficulties in formation of plan- 
etesimals from dust grains,'^ discussion here starts with planetesimals, assuming 
that they have been formed by some mechanism. Planetesimals are bound by their 
self-gravity, but not by material strength. Except for the last stages of planet for- 
mation in which large planets perturb planetesimals, the relative velocities between 
planetesimals do not exceed the escape velocity from their surface, so collisions be- 
tween them mostly result in accretion rather than bouncing or disruption. Since 
orbital damping timescale due to drag forces from disk gas and mean collision/close 
encounter timescale are much longer than orbital periods, planetesimal motion is 
well described by the Kepler motion in most of time. The Keplerian orbital elements 
such as semimajor axis, orbital eccentricity, and orbital inclination are invariable in a 
two-body problem. These elements are changed impulsively by occasional collisions 
and close encounters and gradually by gas drag, although secular changes due to 
secular perturbations are sometimes important. 

Mass distribution of planetesimals evolves due to coagulation. While the col- 
lisional cross-section is determined by the relative velocity, the relative velocity is 
regulated by the mass distribution. This complicated interplay leads to the non-linear 
planetesimal growth. Nevertheless, these processes are better understood than plan- 
etesimal formation from dust grains and gravitational interactions between a planet 
and a gas disk. 

Since the sun is one of G dwarf stars that commonly exist in the Galaxy and 
planets are considered to be by-products of star formation, it was expected that plan- 
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etary systems commonly exist around other stars in the Galaxy. Although searches 
for extrasolar planets (planets orbiting other stars) started in 1940's, they ended 
without any detection or with false detection until the end of the last century. The 
first extrasolar planet around a main-sequence star, 51 Peg b, was eventually discov- 
ered in 1995, by radial velocity survey detecting wobble motion of the host star due to 
orbital motion of a planet.!^ The discovered planet was a gas giant of M ~ lO^M^, 
however, its orbital radius was as small as 0.05 AU and rotates around the host star 
in only 4 days. Such close-in giant planets are called "hot jupiters." As of 2012, 
more than 750 extrasolar planets have been discovered. In 2011, more than 2000 
candidates have been identified by Kepler space telescope that detects eclipses of 
host stars by planets (transit survey) .13 

Most of the extrasolar planetary systems discovered so far have architectures 
quite different from the solar system.--' Other than hot jupiters, many gas giants 
in eccentric orbits have been discovered. The origin of the diversity of planetary 
systems is now actively investigated.'^ Both the radial velocity and transit sur- 
veys detect planets with relatively small orbital radii. Even with such restriction, 
the current surveys suggest that more than 20% of solar-type stars harbor close-in 
rocky/icy planets of M < 30M^. The ubiquity of rocky /icy planets is another im- 
portant subject. It is inferred that rocky/icy planets commonly exist in habitable 
zones that are the ranges of orbital radii where liquid water can exist on planetary 
surfaces. Thus, possible detection of extrasolar life, in particular, biomarkers on 
the extrasolar habitable planets through astronomical observations are also actively 
being discussed.'^ 

Here we mainly describe the basic dynamical and accretionary processes of plan- 
etesimals in the standard scenario of planet formation by showing A'^-body simula- 
tions. Planetesimal accretion that is regulated by planetesimal dynamics is an impor- 
tant stage of planet formation since it determines the timescale of planet formation 
and the basic architecture of the planetary system. We also briefly summarize the 
dynamical effects of disk gas on planets and the core accretion model for formation of 
gas giants and discuss the diversity of planetary systems, extending and generalizing 
the model to extrasolar planetary systems. The model of protoplanetary disks is in- 
troduced in section [2j We describe planetesimal dynamics and accretion in sections 
[3] and [H The formation of terrestrial planets is presented in section [5l We introduce 
the current understanding of the planet-gas disk interaction and the core accretion 
model in section O Section [7] is devoted for a summary. 
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§2. Protoplanetary Disk Model 

To investigate planet formation from a protoplanetary disk, we need a model 
of protoplanetary disks. The standard disk model for solar system formation is the 
minimum-mass solar nebula (MMSN) model that is inferred from the mass distribu- 
tion of planets in the present solar system in which the disk surface density (vertically 
integrated density) scales with the heliocentric distance as r~^/^P 

Unfortunately current radio observations do not have enough resolutions to de- 
termine the surface density distribution of protoplanetary disks although a new large 
radio interferometer in Chile, ALMA, will bring us a lot of information on proto- 
planetary disks (e.g., Refs. [Tn|) . [TT]) ). Thus, we adopt a generalized disk model that 
is based on the MMSN model. The surface density distributions of dust (rock/ice) 
and gas components are given by 

/ r \-3/2 „ 
i:d = 10eice/d -TFT gem ' (2-1) 



lAU. 



i;g = 2400/g(^^j gcm-^ (2-2) 

where /d and /g are multiplicative factors to scale disk surface densities of dust and 
gas components. The disk with /g = /d = 1 corresponds to a 50% more massive 
MMSN model. 

The step function eice is 1 inside the ice line Oice and 4.2 outside a^^e- The 
ice line is the location where the disk temperature coincides with the condensation 
temperature of H2O. In this simple prescription, we set the ice line to that determined 
by an equilibrium temperature in optically thin disk region^ 

where and Lq are the stellar and solar luminosity, as 

/r \i/2 

aice = 2.7 AU. (2-4) 

Note that the magnitude of eicc may be modified by the local viscous dissipation 
and stellar irradiation. Although i7d and are not necessarily proportional to r~'^/^ 
in general, we here present results with this dependence for simplicity (for different 
r-dependence models, see Ref|12j)). The total mass of MMSN within radius ~ 30 
AU is ~ O.OIM0 and radio observations show that the total disk mass around young 
stars is distributed in a range of 0.001-0. lM0.[i3 So, the range of /d and /g may 
be 0.1-10. Hereafter we call the disk model with /g = /d = 1 as the standard disk 
model. 

The rotational angular velocity of gas taking into account the radial pressure 
gradient dP/dr of disk gas is 



^^=^^ + -57) -^(1-^), (2-5) 
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where 

Cg is the sound velocity of gas, f2 and vk are the Kepler angular velocity and velocity, 
and we used P = c^p!^ 
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§3. Planetesimal Dynamics 

The timescale of planetesimal accretion is usually much longer than that of 
planetesimal dynamics in which planetesimal orbits evolve by mutual gravitational 
interaction. Thus planetesimal dynamics controls planetesimal accretion. In this 
section we review the basic planetesimal dynamics due to two-body gravitational re- 
laxation. The two basic effects of two-body gravitational relaxation in a planetesimal 
disk are viscous stirring and dynamical friction: viscous stirring increases the random 
velocity of planetesimals, while dynamical friction realizes the energy equiparation 
of the random energy. We also introduce the orbital repulsion of protoplanets. 

3.1. Two-body relaxation 

The equation of motion for a planetesimal is given as 

^ = -™*^V+ ^ G-.^r5l + ^.as + Fcol, (3-1) 

where m and x are the mass and position of planetesimals, Af^, is the stellar mass, 
and G is the gravitational constant. The terms on the r.h.s. express from left to 
right the solar gravity, the mutual gravitational interaction of planetesimals, the 
force from disk gas, and the velocity change by collisions. We neglect the indirect 
term due to planetesimals since the total mass of planetesimals is much smaller than 
the solar mass. For a planetesimal disk, the solar gravity is dominant except for the 
rare cases of close encounters of planetesimals and thus the orbit of planetesimals is 
almost Keplerian. The mutual gravitational interaction of planetesimals is the main 
perturbing force in a planetesimal disk. The effect of gas will be discussed in section 

El 

The orbit of planetesimals is characterized by the semimajor axis a, the eccentric- 
ity e, and the inclination i. The orbital eccentricity and inclination of planetesimals 
increase as a result of gravitational scattering among planetesimals. The deviation 
velocity, v, of a planetesimal from the Kepler velocity vk of the local non-inclined 
circular orbit is called as random velocity, which is given by 

i; ~ (e^ +i2)^/2fK. (3-2) 

The random velocity is an important factor that controls planet formation as will 
be shown later. For example, the growth timescale of planetesimals depends on the 
random velocity as tgrow oc v'^ when gravitational focusing is effective in collisions. 
Note that the thickness of a planetesimal disk is proportional to the inclination 
(random velocity in a vertical direction) of planetesimals. 

The orbit of planetesimals gradually changes due to the mutual gravitational 
interaction. This process is equivalent to the two-body relaxation process in star 
clusters. In terms of stellar dynamics, a planetesimal disk is a collisional system like 
globular clusters, in the sense that the system evolves through two-body encounters. 

The timescale of two-body relaxation for an equal-mass (m) many-body system 
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is given by 

W = - ' = (3-3) 

where n is the number density of constituent particles, a is the velocity dispersion, 
Tg is the gravitational radius given by Vg = Gm/a"^, and InTl is the Coulomb log- 
arithm.131' The denominator WKr^a means the number of close encounters with 
deflection angle > 90° per unit time and In A takes into account the effect of distant 
encounters. On this timescale a particle forgets the initial orbit. 

The two roles of two-body relaxation in a planetesimal disk are viscous stirring 
and dynamical friction. In the following, we illustrate these processes by showing ex- 
amples of A^-body simulations. The initial conditions of planetesimals used below are 
summarized as follows: 1000 equal-mass (m = 10^^ g) planetesimals are distributed 
in a ring of the radius a = 1 AU with width Aa = 0.07 AU. The surface density 
of the ring is consistent with the standard disk model. The initial distributions of 
eccentricities and inclinations are set by the Rayleigh distributions with dispersions 
(e2)V2 = ^^2^1/2 ^ 2rH/a, where rn is the Hill (Roche) radius of planetesimals given 
by 

The Hill radius is the radius of the gravitational potential well of a body in a rotating 
frame. Note that (e^)^/^ and (i^)^/^ are proportional to the velocity dispersions 
in radial and vertical directions, respectively, and o ~ ((e^) + (i^))"'^/^fK- We set 
= Mq. 

The basic gravitational relaxation processes such as gravitational scattering 
among planetesimals are scaled by the Hill units (rn and 17"^).'^^ 

3.2. Viscous stirring 

The mutual gravitational interaction increases the random velocity of planetesi- 
mals on average. In other words, the Kepler shear velocity due to differential rotation 
is transfered to the random velocity. This process is called as viscous stirring. 

Here we consider an equal-mass planetesimal disk for simplicity. In the dispersion- 
dominated regime ((e^)^/^, (i^)^/^ > rn/a) where the relative velocity of planetesi- 
mals is mainly determined by the random velocity, the viscous stirring rate for (e^) 
is given by -*^^ 

~ ifl . (3.5) 

dt ^relax 

We can obtain the viscous stirring rate for (i^) by replacing e with i in ()3-5p . 

Figures [1] show the system snapshots at t = and 10000 year on the a-e and 
a-i planes. The eccentricities and inclinations of most planetesimals significantly 
increase in 10000 years. On average the increase of e is larger than that of i. The 
distributions of e and i relax into the Rayleigh distributions. We also see the diffusion 
of planetesimals in a, which is the result of random walk in a due to two-body 
scattering. These increases of e and i with the diffusion in a are the basics of viscous 
stirring. 
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Fig. 1. Snapshots of the planetesimal system on the a-e (left) and a-i (right) planes at t = (top) 
and 10000 year (bottom). 




Time evolutions of (e^)^/^ and (i^)^/^ are shown in Figure [2j It is clearly shown 
that (e^)^/^ and (i^)^/^ increase with time as t^/^, and {e^y/'^ / (i^)'^/'^ ~ 2 that 
corresponds to the anisotropic velocity dispersions in radial and vertical directions. 
These properties of {e^Y^"^ and (i^)-*^/^ are the characteristics of two-body relaxation 
in a disk."*^ 

As the number density of planetesimals is inversely proportional to the disk 
thickness that is proportional to the velocity dispersion, n oc a~^, we have treiax oc cj^. 
Thus, we have a oc t^^^ from ()3-3p . The origin of the anisotropy of the velocity 
dispersion is the shear velocity between planetesimals with different a due to the 
differential rotation in the Kepler potential. This type of anisotropy is also known 
for the Galactic stellar disk.'I^ 

In reality the gas drag and collisions damp the random velocity of planetesi- 
mals. In planetesimal accretion, there exist equilibrium values of the eccentricity 
and inclination, at which viscous stirring and damping by the gas drag and collisions 
balance. 
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Fig. 3. Snapshots of the planetesimal system on the a-e (left) and a-i (right) planes at t = (top) 
and 3000 year (bottom). The large circle indicates the protoplanet. 
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Fig. 4. Time evolution of cm (solid) and im (dotted) of the protoplanet. 



3.3. Dynamical friction 

Next we consider the two-body relaxation for particles with a mass distribu- 
tion. Dynamical friction is the process of the equiparation of the random energy, 
{l/2)mv'^ , of planetesimals. In other words, the random velocity becomes v oc m^l"^ . 
As an illustration of dynamical friction, we consider a simple case with a protoplanet 
(large planetesimal) with M = 100m embedded in a swarm of planetesimals. We 
focus on the orbital evolution of the protoplanet. The initial orbital elements of the 
protoplanet are qm = 1 AU and cm = iu = 0.01. 

For M ^ m and cm > e, the dynamical friction rate for cm is approximated 



dp2 



m e 



M 



(3-6) 

dt Mt,el^^ 

The rate for zm is obtained by replacing cm with im in ()3-6p . Figures [3] show the 
system snapshots at t = and 3000 year on the a-e and a-i planes. We see that the 
eccentricity and inclination of the protoplanet decrease to almost in 3000 years. 
However, its semimajor axis is kept almost constant. On the other hand, the eccen- 
tricities and inclinations of the neighbor planetesimals of the protoplanet are raised 
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by reaction. The V-like structure around the protoplanet on the a-e plane corre- 
sponds to the constant Jacobi energy curve. This heating of neighbor planetesimals 
by a protoplanet leads to the decrease of the growth rate of the protoplanet as will 
be shown in section [4.31 ^^ 

Time evolution of e m and i M of the protoplanet is shown in Figure HI In 3000 
years, the eccentricity and inclination of the protoplanet are reduced to ~ 0.001 
due to dynamical friction from small planetesimals, in other words, the orbit of the 
protoplanet becomes a non-inclined nearly circular orbit. 

One of the important features of dynamical friction is that the dynamical friction 
rate does not depend on the mass of individual particles but the system density. For 
protoplanet-planetesimal scattering, trciax for the protoplanet is treiax oc (nAf^)~^ 
from ([331). Therefore, ([Ml) leads to |de|^/dt| oc nm = p, where p is the system 
density. 

When the mass of a protoplanet is smaller than the total mass of planetesimals 
in its feeding zone, the recoil of planetesimal scattering damps cm and iu of the 
protoplanet as explained above. The recoil can also change the semimajor axis gm 
of the protoplanet under some circumstances. We comment on it in section [6l after 
explaining orbital migration of a protoplanet due to interactions with a gas disk. 

3.4. Orbital repulsion 
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Fig. 5. Snapshots of the planetesimal system on the a-e (left) and a-i (right) planes at t = (top) 
and 2000 year (bottom). The large circles indicate the protoplanets. 



As an application of dynamical friction let us think the orbital evolution of two 
interacting protoplanets embedded in a swarm of planetesimals. It is found that if the 
orbital separation, b, of the two protoplanets is smaller than a few times their mutual 
Hill radius, they expand their orbital separation to 6 > Srn-'^ This phenomenon is 
called as orbital repulsion. In Figures [5l the protoplanets with M = 100m initially 
on non-inclined circular orbits with the orbital separation of 6 = 3rH expands the 
orbital separation to 6 ~ Srjj keeping nearly non-inclined circular orbits in 2000 
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Fig. 6. Schematic illustration of orbital repulsion on the a-e plane. 



years. 



The orbital repulsion is a coupling effect of gravitational scattering between large 
bodies and dynamical friction from small bodies. Figure [6] schematically explains the 
mechanism of orbital repulsion. The mechanism consists of two stages: 1) Scattering 
between two protoplanets on nearly circular non-inclined orbits increases their ec- 
centricities and orbital separation. 2) Dynamical friction from planetesimals reduces 
the eccentricities of the protoplanets while the orbital separation remains almost 
the same. Thus, the orbital separation of two protoplanets increases keeping nearly 
circular orbits. If the orbital separation is less than 5rH, relatively strong scattering 
occurs and the separation expands rapidly. We can analytically show the behavior 
of the first stage, based on the conservation of energy and angular momentum in 
two-body scattering under the solar gravity. 

During the course of planetesimal accretion, as protoplanets grow, their orbital 
separation normalized by the Hill radius decreases, since rn oc M^/^. This implies 
that they repeat the orbital repulsion while growing. Consequently, the orbital 
separation of protoplanets is kept always larger than about Srjj. This is one of the 
important factors that realize oligarchic growth of protoplanets. The typical orbital 
separation is 6 ~ lOrn, which only weakly depends on the mass and the semimajor 
axis of protoplanets and the disk surface density.f^S 
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§4. Planetesimal Accretion 

Orbiting the sun planetesimals sometimes collide with each other to form pro- 
toplanets. In this section we review the basic physics of accretionary evolution of 
planetesimals. Runaway growth of planetesimals and oligarchic growth of protoplan- 
ets are demonstrated by showing A^-body simulations. 

4.1. Growth modes of planetesimals 

In general, there are two major modes when particles grow by coagulation, 
namely, "orderly" and "runaway" growth. In the orderly growth mode, all the par- 
ticles grow equally, in other words, mass ratios between particles tend to be unity. 
On the other hand, in the runaway growth mode, larger particles grow faster than 
smaller ones and their mass ratios increase monotonically. Which growth mode is 
relevant to planetesimal accretion had been controversial around the end of the last 
century. 

Now we consider the evolution of mass ratio between two test particles with 
mass Ml and M2 (Mi > M2). The time derivative of the mass ratio is given by 

d {Mi\ Ml ( 1 dMi 1 dM2\ ^^^^ 



dt \M2) M2 \Mi dt M2 dt 

which shows that it is the relative growth rate (l/M)(dM/dt) that determines the 
growth mode of particles. If the relative growth rate decreases with M, d{Mi/M2) /dt 
is negative and the mass ratio tends to be unity. In this case, the growth mode is or- 
derly. On the other hand, if the relative growth rate increases with M, d{Mi/M2) /dt 
is positive and the mass ratio increases, which results in runaway growth. 

The growth rate of a test planetesimal with mass M and radius R by accreting 
field planetesimals with mass m (M > m) is given by 

,2 



^ ~ UmTTR^ (1 + ^] v,,,m, (4-2) 

where is the number density of the field planetesimals, and frei and Vesc are the 
relative velocity between the test and the field planetesimals and the surface escape 
velocity from the test planetesimal, respectively.!^ For simplicity, we assume that all 
collisions lead to accretion (sticking probability is unity). Here 7ri?^(l + fesc/^rei) 
the collisional cross-section with gravitational focusing, which can be easily obtained 
using the conservation of energy and angular momentum. The term I'esc/^rcl indicates 
the enhancement of geometrical collisional cross-section by gravitational focusing. 
We discuss the growth mode of planetesimals using ()4-2p below. 



4.2. Runaway growth of planetesimals 

In the early stages of planetesimal accretion, the growth mode of planetesimals 
is runaway growth, where larger planetesimals grow more rapidly than smaller ones 
and their mass ratios increase with time.'^''^''^ We illustrate the basic processes 
of runaway growth of planetesimals by showing the results of A^-body simulation 
of planetesimal accretion. Figures [7] and [5] show an example of runaway growth.^^J 
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Fig. 7. Snapshots of the planetesimal system on the a-e plane for t = 0, 10^, and 2 x 10^ year. 
The circles represent planetesimals and their radii are proportional to the radii of planetesimals. 



Planetesimals are initially distributed in a ring around a = 1 AU with the width 
Aa = 0.02 AU. The planetesimal system initially consists of 3000 equal-mass (m = 
10^^ g) bodies, which is consistent with the standard disk model. The internal density 
of planetesimals is /Op = 2 gcm~^. For simplicity, the calculation in this section is 
gas- free. 

Figure [7] shows the system snapshots on the a-e plane for t = 0, 10^, and 2 x 10^ 
year. In 2 x 10^ years, the number of bodies decreases to 1322. It is clearly shown 
that a planetesimal grows in the runaway mode. At t = 2 x 10^ year, the mass 
of the largest body reaches about 200 times the initial mass, while the mean mass 
of planetesimals becomes only about twice larger. Note that the eccentricity and 
the inclination of the largest body is always kept small due to dynamical friction 
from smaller bodies. These small eccentricity and inclination facilitate runaway 
growth.E3.E3 

The evolution of the mass distribution of planetesimals is shown in Figure [8l 
The cumulative number of bodies Uc against mass is plotted. The mass distribution 
relaxes to the distribution that is well approximated by a power-law distribution. 
The runaway body at t = 2 x 10^ year is shown by a dot that is separated from the 
continuous mass distribution. The mass range lO^^g ^ m < lO^^g, which contains 
most of the system mass, can be approximated by dnc/dm oc m". The power index 
calculated by using the least-square-fit method is a ~ —2.5. This index can be 
derived analytically as a stationary distribution.ESJ The power index smaller than 
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Fig. 8. Snapshots of the cumulative mass distribution of the planetesimal system for t = 
(dotted), 10^ (dashed), and 2 x 10''' (sohd) year. The dot stands for a runaway body at t = 2 x 10^ 
year. 



-2 is a characteristic of runaway growth, which means most of the system mass 
exists in small bodies. Note that runaway growth does not necessarily mean that 
the growth time decreases with mass, but does mean that the mass ratio of any 
two bodies increases with time. The runaway body keeps growing and then isolates 
from the continuous power-law mass distribution. In this stage, the runaway body 
predominantly grows as a sink of the mass flow from the continuous power-law mass 
distribution. 

The runaway growth of planetesimals is explained as follows. When gravitational 
focusing is effective, ^esc/'^rei ^ (|4"2p reduces to 

dM 



dt 



oc 



U^M^I'^v-^, (4-3) 



where i7d is the surface density of planetesimals and we used rim oc Z'df"^, v^sc oc 
Af -"^Z^, R oc M"*^/^, and Wrei — v- On the early stage of planetesimal accretion, and 
V barely depend on M, in other words, the reaction of growth on and v can be 
neglected since the mass in small planetesimals dominate the system. In this case 
we have 

~ c M'/3, (4,4) 

M dt ^ ' 

which leads to runaway growth as shown in section 14. li 

4.3. Oligarchic growth of protoplanets 

Protoplanets are formed through runaway growth of planetesimals. In the late 
runaway stage, protoplanets grow while interacting with one another. Kokubo and 
Ida'^ investigated this stage and found oligarchic growth of protoplanets: similar- 
sized protoplanets grow keeping their orbital separation larger than about Srn, 
while most planetesimals remain small. Through oligarchic growth, a bi-modal 
protoplanet-planetesimal system is formed at the post-runaway stage. 
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Fig. 9. Snapshots of the planetesimal system on the a-e plane for t — 0, 10^, 2 x 10^, and 
4 X 10^ year. The circles represent planetesimals and their radii are proportional to the radii of 
planetesimals. 



We present an example of the A^-body simulation that shows oligarchic growth.l^ 
Figure shows the result of a large scale simulation starting from = 10000 equal- 
mass (m = 1.5 X lO^^g) planetesimals in 0.5 AU < a < 1.5 AU, which is consistent 
with the standard disk model. In this calculation, the 6-fold radius of planetesimals 
is used to accelerate the accretion process. The use of the 6-fold radius of plan- 
etesimals does not change the growth mode of planetesimals but shorten the growth 
timescale about 6 times .^^D 

The accretion propagates from small to large a. This is because the accretion 
timescale is smaller for smaller a since the surface number density of planetesimals is 
higher and the orbital period is smaller for smaller a [see (|4-8|) ]. In 4 x 10^ years, the 
number of bodies decreases to 333. About 10 protoplanets form with mass ~ lO'^^ 
g on nearly circular non-inclined orbits with the orbital separation of ~ lOrn- Note 
that at large a the protoplanets is still growing. The result of the A^-body simulation 
is consistent with the estimation based on the oligarchic growth model described 
below. 

Oligarchic growth is the result of the self-limiting nature of runaway growth of 
protoplanets and orbital repulsion among protoplanets. The formation of similar- 
sized protoplanets is explained by the slow-down of runaway growth.E3'[i9 Once the 
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mass of the runaway body isolates from the continuous power-law mass distribution 
of planetesimals (see Figure [8]) , we can approximate the system by a two-component 
system: a protoplanet and small planetesimals. When the mass of the protoplanet M 
exceeds about 100 times the typical mass of planetesimals, the protoplanet effectively 
heats up the random velocity of neighbor planetesimals to be t; oc 

MV3.1I3 In this 

case, from (|4-3p . the relative growth rate becomes 

' oc E,M-'/\ (4-5) 



M dt 



On this stage, decreases through accretion of planetesimals by the protoplanet 
as M increases!^ That is, the relative growth rate is a decrease function of M. 
Thus the growth mode among protoplanets becomes orderly, in other words, the 
mass ratios between protoplanets tend toward unity rather than increase. Thus, the 
neighboring protoplanets grow keeping similar masses. 

Note that still in this stage, the mass ratio of a protoplanet to its neighbor plan- 
etesimals increases since for its neighbor planetesimals with mass m, (l/m)dm/dt oc 
i7dm-'^/^M~^/^, which leads to 

(l/M)(dM/dt) (mY' 

(l/m)(dm/dt) \m) ' ^ ' 

The relative growth rate of the protoplanet is by a factor of {M/m)^^^ larger than 
that of the planetesimals. Thus, a bi- modal protoplanet-planetesimal system is 
formed. 

While the protoplanets grow, orbital repulsioiJ^ keeps their orbital separations 
~ lOrn as shown in section [37^ 



4.4. Isolation mass of protoplanets 

As a result of the oligarchic growth of protoplanets, a protoplanet system is 
formed from a planetesimal disk. The dynamical characteristics of the protoplanet 
system are estimated based on the oligarchic growth model. In oligarchic growth, 
protoplanets are formed with a certain orbital separation. Given this orbital sepa- 
ration 6, the isolation (final) mass of a protoplanet at a is estimated bj^^J 

M,„ . 2™K. = OMfr^^: ( {^f" ( ^V'^Me, (4-7) 



lOrny VlAU/ V^© 
and the growth timescale of the protoplanet is 
M 

^grow 



dM/dt 

years. (4-oj 



lOiSg / V M, 



Dynamics and Accretion of Planetesimals 



17 



Then the growth timescale of the protoplanet to the isolation mass is given by 

3/5 



W(^iso) ^ 0.7 X loVd '^'/g~'/'e,^e^' 



m 



2/15 
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2gcm 
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b V^^° ( a \ 59/20 
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Fig. 10. Isolation mass of protoplanets against the semimajor axis with the mass of planets in the 
solar system and the ice line (dotted). 

Figure [10] shows the isolation mass of protoplanets against the semimajor axis 
for the standard disk model for solar system formation. In the terrestrial planet 
region Miso ~ O.IM^, which is as large as the mass of Mercury and Mars. This 
suggests that they are leftover protoplanets. In order to complete Venus and Earth 
whose masses are one order of magnitude larger than that of protoplanets, further 
accretion of protoplanets is necessary. In the standard scenario of terrestrial planet 
formation the final stage is the "giant impact" stage of protoplanets. The isolation 
mass of protoplanets in the gas giant region is around S-IOM®, which may accrete 
gas from the disk before the dispersal of disk gas. The mass (~ lOM®) and orbital 
separation (~ 10 AU) of ice giants are consistent with the oligarchic growth model. 
However, their growth timescale is longer than the age of the solar system. 

The oligarchic growth model of protoplanets is now generally accepted as the 
standard process of planet formation though it still has some discrepancies. The gen- 
eralized oligarchic growth model is used to study the diversity of extrasolar planets 
together with the core accretion model as will be shown in section 16.21 
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§5. Terrestrial Planet Formation 

It is generally accepted that the final stage of terrestrial planet formation is 
the giant impact stage where protoplanets (planetary embryos) formed by oligarchic 
growth collide with one another to complete planets;^-''-- This stage is being ac- 
tively studied as many small extrasolar planets are discovered. 

5.1. Orbital instability of protoplanet systems 

10' 

io« 

10» 

•-^ 10* 

m 
C 

-o 

1000 
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10 

2 4 6 8 10 

^ini(^H) 

Fig. 11. Timescale of orbital instability against the initial orbital separation of protoplanets. The 
circles, triangles, squares and crosses indicate the initial (e^)^^^ = 2(i^)^''^ = 0, lh,2h, and 4/i, 
where h is the reduced Hill radius r-a/a. The solid line shows the result of Chambers et al.^^ 

Through oligarchic growth protoplanets with orbital separation ~ lOrn are 
formed from planetesimals. Though protoplanets perturb each other, the proto- 
planet system is orbitally stable when disk gas exists since its gravitational drag 
damps their eccentricities (see section 16. ip . However observationally it is inferred 
that disk gas depletes on the timescale of 1-10 million years.'^ Thus in the long-term 
the protoplanet system becomes unstable through mutual gravitational perturbation 
after the dispersal of the gas disk. 

The timescale of the orbital instability of a protoplanet system is numerically 
obtained by A^-body simulations as 

logtinst ^ Ci ( — ] +C2, (5-1) 

where 6ini is the initial orbital separation of adjacent protoplanets and ci and C2 
are functions of the initial (e^)^/^ and (i^)^/^ of the system.l^'IS^ The physical 
interpretation of this semi-logarithm dependence is still unclear. 

Figure [TT] shows the orbital instability timescale of a protoplanet system consist- 
ing 10 equal-mass (O.IM^) protoplanets obtained by A^-body simulations.—J Here 
"instability" means orbital crossing or collision of protoplanets. The semi-logarithm 
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dependence on the initial orbital separation of protoplanets is clearly shown. The 
constants ci and C2 are smaller and larger for the larger initial (e^)^/^ and 2(i^)^/^ 
of protoplanets, respectively. With only a small initial (e^)^/^ and 2(i^)^/^ the in- 
stability timescale drastically shortens for large 6ini. 

5.2. Giant impacts 




2x10' 4x10' 6x10' 8x10' 10» 0.5 1 1.5 8 2.5 



t(yr) a(AU) 

Fig. 12. Time evolution of the semimajor axes (solid lines) and pericenter and apocenter distances 
(dotted lines) of planets (left). Snapshots of the system on the a-e plane a,t t — 0, 10^, 10^, 10*, 
and 2 x 10* year (right). The size of circles is proportional to the physical size of planets. 

After a protoplanet system becomes orbitally unstable the giant impact stage 
of protoplanets begins. As this process is stochastic in nature, in order to clarify 
it, it is necessary to quantify it statistically. Kokubo et al.^ Kokubo and Ida,'32l' 
and Kokubo and GendsP^ investigated the basic dynamics of the giant impact stage 
statistically with many A^-body simulations. 

Figure [12] shows an example where three terrestrial planets are formed from 16 
protoplanets.'^ In this run the giant impact stage lasts for about 7 x 10'' years. 
In the standard disk model, two Earth-sized planets typically form in the terrestrial 
planet region. Kokubo et al.® find that an important parameter of the initial 
protoplanet system for the number and mass of final planets is the total mass of 
protoplanets, Mtot- The effects of the surface density distribution of the disk are 
unified using Mtot ■ In Figure [T31 the average masses of the largest planet (Mi ) 
and the second-largest planet (M2) are plotted against Mtot for various models of 
protoplanet systems with different surface densities and radial extent together with 
their empirical fits. It is clearly shown that both (Mi) and (M2) increase almost 
linearly with Mtot and (Mi) ~ 0.5Mtot and (M2) ~ O.SMtot- This result shows that 
protoplanet accretion proceeds globally, in other words, over the whole terrestrial 
planet region. Thus the large-scale radial mixing of material is expected. 

The spin parameters of terrestrial planets are determined by the angular mo- 
mentum brought by giant impacts-I^SJ The spin angular velocity averaged in mass 
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Fig. 13. Average masses of the largest (Mi) (filled symbols) and second largest {M2) (open 
symbols) against Aftot and the empirical fits for (Mi) (solid line) and (M2) (dotted line). The 
triangles, circles, and squares correspond to /d = 0.3, 1, and 3 models with flout ~ 1.5,2.0,2.5, 
and 3.0 AU. 



bins against planet mass is shown in Figure [Hi.Ef Kokubo and GendeP used 
the realistic accretion condition of protoplanets obtained by the numerical collision 
experiments of protoplanets.'® It is clear that the average angular velocity is al- 
most independent of mass. The average values are about 70% of the critical angular 
velocity for rotational breakup 

U^cr = (^-^J . (5-2) 

This is a natural outcome for the giant impact stage where the impact velocity of 
protoplanets are ~ Vesc- For an Earth-mass planet, the spin angular velocity of 
O.Twcr corresponds to the spin angular momentum of 9.5 x 10^^ g cm^ s~^, which 
is 2.7 times larger than that of the Earth-Moon system. So the angular momentum 
of the Earth-Moon system is not a typical value but it is reasonably within the 
distribution of uj. 

In Figure [Tib, we show the obliquity distribution with the isotropic distribution 

nde = -sinede. (5-3) 

The obliquity ranges from 0° to 180° and follows an isotropic distribution. The 
isotropic distribution of e is a natural outcome of giant impacts. During the giant 
impact stage, the thickness of the protoplanet system is ~ a{i'^)^^'^ ~ lOrn, which is 
much larger than the size of protoplanets R ~ 10~^rH, where i is the inclination of 
protoplanets. Thus, collisions are three-dimensional. We find that these results are 
independent of the initial protoplanet system parameters. 

Our results clearly show that the distribution of spin obliquity is isotropic, which 
suggests that terrestrial planets formed through giant impacts are likely to have 
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Fig. 14. Average spin angular velocity against planet mass (left) and normalized cumulative 
distribution of obliquity with the isotropic distribution (dotted line) (right) of all planets formed 
in the 50 runs of the standard model. The error bars indicate the standard deviation and the 
dotted line shows ujcv (left). 



£ ~ 90°, in other words, their spin axes are nearly on their orbital plane. Prograde 
spin with small obliquity, which is common to terrestrial planets in the solar system 
except for Venus, is not a common feature for planets assembled by giant impacts. It 
should be noted, however, that the initial obliquity of a planet determined by giant 
impacts can be modified substantially by stellar tide if the planet is close to the star 
(Mercury) and by satellite tide if the planet has a large satellite (Earth) . 
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§6. Planet-Gas Disk Interaction 

So far we have discussed dynamics and accretion of solid (rocky and icy) planets. 
In the solar system, gas giant planets, Jupiter and Saturn, exist. Since the gas giants 
are 100 times or more massive than terrestrial planets, the architecture of planetary 
systems is sculpted by the giants, if they exist in the systems. The mass of protoplan- 
etary disks is about one order of magnitude larger than that of giant planets, which 
implies that planetary orbits are significantly affected by interactions with disk gas. 
Here we briefly summarize current understanding of planet-gas disk interactions and 
formation process of gas giant planets. Since distributions of gas giants in extrasolar 
planetary systems have been revealed by rapidly developing observations, we can 
compare the theoretical predictions with the observed distributions. Although some 
aspects of the observed distributions are explained by the theory, there still remain 
important unsolved problems. 

6.1. Orbital evolution 

While orbits of planetesimals and protoplanets can be eccentric, disk gas rotates 
in a circular orbit. Thereby, planet-disk interaction generally damps the orbital 
eccentricity of planetesimals and protoplanets. For small planetesimals, aerodynamic 
gas drag is dominant because the ratio of the surface area to the volume is higher for 
smaller bodies. On the other hand, dynamical friction from disk gas is dominant for 
large bodies (protoplanets) that are more massive than the lunar mass O.OIM®). 

The damping timescale of the relative velocity u between a planetesimal (mass 
M and physical radius R) and disk gas due to aerodynamical gas drag, in the case 
where Reynolds number is larger than unity, is given 

Mu 

The relative motion between a planetesimal/protoplanet and disk gas has two com- 
ponents: a component due to an eccentric orbit of the planetesimal (~ cvk) and an 
offset due to the pressure gradient in the gas (~ tjvk)- Then, we have u ~ (e + r])vK- 
On a timescale of tstop, both evK and tjvk are damped. Thus, eccentricity 
damping timescale is 

_ e 

where the gas density pg is evaluated by Ug divided by disk thickness ~ 2cs/f2, and 
the sound velocity Cg is given by the equilibrium disk temperature T. 

The damping of tjvk causes inward radial migration of a planetesimal, because 
the planetesimal is subject to "head wind" to loose angular momentum, L. Thus, 
dL/dt ~ -Mar]VK/tstop- Since dL/dt ~ Md{avK)/dt = (l/2)Mt>K(da/dt), the 



(6-1) 
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migration timescale is given by 

— ^ 

Equations (|6-2p and (j6-3p show that tdrag.a is 2-3 orders of magnitude longer than 

^drag,e- 

The eccentricity damping timescale due to dynamical friction from disk gas is 
given by 

tdi ~ — , (6-4) 

where re = GM/c^ is the Bondi radius of the protoplanet. The last factor in the 
denominator is Cg, which means that the velocity of density wave propagation is ~ Cg 
rather than n, since we consider the case of Cs > u. Note that the final form of tdi 
is similar to trelax (|3'3p if a is replaced by Cg. The eccentricity damping timescale is 
given bjJSSJ 

, 1 / M \ r \3/2 /M* \3/2 

~ 250/-M— (— -) -1 years. (6-5) 



^M©y VlAU/ V^^0, 

Equations (|6-ip and (|6-4p show that the effects of gas drag and dynamical fric- 
tion are weaker for larger and smaller bodies, respectively, and e-damping is the 
weakest for M ~ O.OIM^ at r ~ 1 AU. The e-damping timescale even for bodies 
with M ~ O.OIM^ is shorter than the disk lifetime ~ 1 — 10 million years. While 
the eccentricities of small planetesimals are excited by nearby protoplanets, allowing 
mutual collisions and accretion by the protoplanets, the eccentricities of the pro- 
toplanets are not excited, so that protoplanet-protoplanet collisions are generally 
inhibited in the presence of gas. 

Although planet-disk gravitational interactions are much more complicated than 
aerodynamical gas drag, the timescale of radial migration due to dynamical friction 
is evaluated in a similar way as 

1 , , 

td{,a ~ 7^*df- (6-6) 
zr] 

This migration is often called "type-I migration." Tanaka et al.'^ derived a numer- 
ical factor for ()6-6p through detailed 3-D linear calculation, taking into a curvature 
effect (geometrical imbalance between the effects from inner and outer disks). With 
the numerical factor, (16-611 is 
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The migration timescales for an Earth-mass planet at lAU and a core of ~ lOM^ 
at 5 AU are ~ 10^ years, which are 10-100 times shorter than disk lifetime, sug- 
gesting that the Earth and Jupiter's core cannot survive. This is one of the most 
serious difficulties in planet formation theory today. On the other hand, discov- 
ery of many close-in planets in extrasolar planetary systems requires this kind of 
inward migration. Many possibilities to retard, halt or reverse the migration are 
discussed, although the mechanism for planets to survive type I migration has not 
been identified yet. 

Orbital migration of protoplanets due to planetesimal scattering (gravitational 
interactions with a planetesimal disk) also exists, as pointed out in section 13. 3i 
Although there is no systematic effect like r] in the gas disk, a curvature effect 
still exists. While i7d <C Sg, Cs is generally larger than the velocity dispersion of 
planetesimals, so that "planetesimal-driven migration" cannot be always neglected 
compared with type I migration.'SS' However, the migration speed and direction have 
not been clarified yet. 

6.2. Formation of gas giants 

For the formation of gas giant planets, two models have been proposed. In the 
"core accretion" model, a rocky/icy core first accretes from planetesimals and gas 
accretion onto the core follows.'^ In the "disk instability" model, on the other 
hand, a protoplanetary disk fragments due to gravitational instability to directly 
form a gas giant (s).!^ Since rocky/icy planets cannot be formed by the disk insta- 
bility and even Jupiter and Saturn are metal-rich compared to the solar composition, 
the core accretion model has been regarded as a "standard" model. Here, we sum- 
marize the core accretion model. Note that a possibility of disk instability is now 
being revisited, because extrasolar gas giants with large orbital radii, which are not 
easy to be formed in situ through the core accretion, are being observed by direct 
imaging. 

When the core mass becomes larger than a critical core mass, pressure gradient 
no more supports envelope gas hydrodynamically against the core's gravity and 
hydrostatic envelope no longer exists. The critical core mass (Mc^cr) depends on 
planetesimal accretion rate onto the core (Mc) and the grain opacity in the envelope 
(Kgr). Through 1-D calculation, Ikoma et al^^ found that 

Mcc^iof -1 (^) (6-8) 

where Kgj-^p (~ icm^g"^) is the grain opacity given by Pollack et al.,'33' who assumed 
dust grains with interstellar abundance and size distributions. Faster accretion and 
higher opacity (relatively large Mc and Kgr) result in a warmer planetary envelope 
and enhanced pressure gradient, so that Mc^cr increases .^^'H^ 

After Mc exceeds Mc^cd heat generation due to gravitational contraction of gas 
envelope itself supports the envelope against dynamical collapse and the envelope 
undergoes quasi-static contraction.HB' The contraction allows disk gas to flow from 
the protoplanetary disk into the Bondi radius of the planet, so that the contraction 
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rate determines the rate of gas accretion onto the planet. 

Although the rate of the quasi-static contraction of envelope should be evaluated 
by radiative transfer calculation, we here present a simple estimate for the rate.!^ 
For Mc ~ Mc^CT) heat generation due to planetesimal accretion marginally equili- 
brates with the gravity of the core. In the quasi-static envelope contraction, heat 
generation due to the envelope contraction marginally equilibrates with the gravity 
of the planet. Replacing and Mc in ()6-8p by M (total mass of the core and 
envelope) and M/ig^acc, the gas accretion timescale tg,acc is given by 



''g,acc 



10^ 



M 



-(2.3-4) 



%r,P 



years. 



(6-9) 



Detailed radiative transfer simulations^ 'SSJ showed consistent results. 

Since tg^acc rapidly decreases with increasing M, the gas accretion is a runaway 
process and initial contraction regulates the total gas accretion timescale. Equation 
(|6-9|) shows that for gas accretion to start the runaway gas accretion within disk 
lifetime (tdep) of 1-10 million years, the core mass more than a few to ten Earth 
masses is required for Kgr ~ (O.Ol-l)Kgr^p. 

6.3. Formation sites of cores 

Using the results in section 14.41 and 16. 2^ we can quantitatively determine the 
formation sites of cores with a few or more, neglecting type-I migration, which 
is highly uncertain. We consider a conservative condition for formation of gas giants 
such that Mc > lOM®. In the following we use b = lOrn, Pp = 2.0 g cm^'^, and 
m = 10^® g. 

In the inner regions, cores acquire the isolation mass Miso given by (|4-7p . which 
shows that Miso > IOM0 only at r > oin where 
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(6-10) 

where we assumed oc M^. This equation shows that preferred locations of cores 
for giant planets are the regions beyond the ice line (r > Oice), although formation 
of large cores is possible even inside the ice line for the most massive disks with 
/d > 10. 

On the other hand, in outer regions, core growth is so slow that the core mass 
does not reach the isolation mass before the depletion of disk gas. Then, the condition 
for formation of cores with Mc larger than lOM^ is tgrowllOM^) < t^ep, where we 



assume tdep ~ 10 million years and tgrow is given by 



This condition imposes 
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an outer limit of the semimajor axis (oout) for formation of gas giants: 
. ,,14/27 /eiceV°/'V*^*A'^'' / tdep V'°^'' ATT 

where solar abundance, /d = /g, is assumed. Formation of cores for gas giants is 
possible in the regions with a\a < r < Oout- 

6.4. Diversity of planetary systems 
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Fig. 15. Schematic illustration of diversity of planetary systems. 

Using the results in the last subsection, we discuss the diversity of planetary 
systems.liS'Ef (Figure US]). We now know that M > IOM0 within tdep only at 
Oin < < Oout- In low-mass disks for which Oin > Cout, gas giants are not formed. 
Equations (|6-lUp and (|6-11|) suggest that formation of gas giants is not expected in 
disks with /d ^ 1 (disks with masses smaller than the mass of the standard disk). On 
the other hand, in high-mass disks, it is expected that multiple gas giants are formed 
since the range between a\a and Oout is broad. The multiple giants often undergo 
close scattering among them, resulting in formation of gas giants in eccentric orbits 
that are found in extrasolar planetary systems. In massive disks, cores and gas giants 
form early enough when full amount of disk gas still exists. Then, the gas giant that 
has opened up a gap in the gas disk would migrate together with disk accretion (type 
II migration j^B'BSJ the vicinity of the host star. Thus, initial disk masses would 
produce diversity of planetary systems.^ 

Gas accretion rate onto the core and the mechanism to stop gas accretion are 
being discussed by many authors, as well as realistic type I migration rate. Although 
these processes sill have large uncertainties, the planetesimal dynamics and accretion 
described here are relatively well understood and they are fundamental pieces for 
discussion on formation of planetary systems. 
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§7. Summary 

In the standard scenario of planet formation terrestrial and icy planets and cores 
of gaseous planets are formed by accretion of planetesimals. The random velocity 
(eccentricity e and inclination i) of planetesimals controls planetesimal accretion and 
then growth of planetesimals changes the random velocity of planetesimals. This 
interplay of planetesimal dynamics and accretion show interesting phenomena of 
planetesimal accretion. We have demonstrated the basic dynamical and accretionary 
processes of planetesimals by showing the examples of A^-body simulations. 

In orbital evolution of planetesimals two-body gravitational relaxation plays key 
roles. The important basic processes are viscous stirring and dynamical friction. 
Viscous stirring increases dispersions (e^)^/^ and (i^)^/^ in proportion to time as 
t^/*^, keeping the ratio {e^)^^"^ / {i'^y^'^ ~ 2, which is the characteristics of two-body 
relaxation in a Keplerian disk. Dynamical friction realizes the equiparation of the 
random energy, in other words, e, i oc m~^/^. Orbital repulsion of two protoplan- 
ets embedded in a swarm of planetesimals keeps the orbital separation of the two 
protoplanets b > Sfh- Orbital repulsion is one of the key processes that realize the 
oligarchic growth of protoplanets. All these elementary processes control the basic 
mode, timescale, and spatial structure of planetesimal accretion. 

On the early stage of planetesimal accretion the growth mode of planetesimals is 
runaway growth where larger planetesimals grow faster than smaller ones. This is due 
to gravitational focusing that enhances the collisional cross-section of planetesimals 
by the self-gravity. The mass distribution of planetesimals is well approximated 
by a power-law distribution of ndm oc m~"dm, where a ~ 2.5 for small bodies, 
and detached runaway bodies. Once the mass of a protoplanet (runaway-growing 
body) exceeds a critical mass, it effectively heats up the random velocity of local 
planetesimals, which results in orderly growth among protoplanets. Then orbital 
repulsion between adjacent protoplanets keeps their orbital separations b ~ lOrn- 
This growth mode is called as oligarchic growth of protoplanets where similar-sized 
protoplanets predominantly grow with certain orbital separations. These processes 
of planetesimal dynamics and accretion are relatively well understood and are now 
incorporated into the standard model of planet formation. 

The final stage of terrestrial planet formation is known as giant impact stage. 
A protoplanet system becomes orbitally unstable after dispersal of disk gas on the 
timescale logtmst — ciftini/^H + C2 where 6ini is the initial orbital separation of pro- 
toplanets and ci and C2 are functions of the initial (e^)^/^ and (i^)^/^ and then the 
giant impacts of protoplanets start. On this stage accretion proceeds globally and 
the total mass of protoplanets Mtot is a key parameter that determines the mass of 
planets. The masses of the largest and the second-largest planets increase with Mtot 
almost linearly. The spin parameters of planets are determined by giant impacts. 
The RMS spin angular velocity oj is as large as 70% of the critical spin angular 
velocity for rotational instability. The obliquity of the planets follow an isotropic 
distribution. 

We also discussed the orbital evolution of planets by the planet-gas disk inter- 
action. The core accretion model for gas giant formation is summarized and used to 
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discuss the diversity of planetary systems. 

As we have shown, planet formation consists of multi-scale, multi-layer pro- 
cesses regulated by a variety of physics. Some of fundamental processes are still 
unclear, in particular, formation of planetesimals from dust and orbital migration 
of (proto)planets. Both processes are regulated by structure and evolution of a pro- 
toplanetary disk. We hope ALMA will reveal detailed structure/evolution of proto- 
planetary disks. The diversity of observed extrasolar planets suggests that secular 
gravitational perturbations among planets also play an important role in creating 
the architecture of planetary systems. We need to deeply explore secular orbital 
dynamics as well as planet-disk interactions to understand the diversity of planetary 
systems. We can calibrate theoretical models by using the data of rapidly developing 
observations of extrasolar planets. In the theoretical model of planet formation, there 
are many unsolved problems related with many different kinds of physics. However, 
we already have tools to attack the problems. Studies on planet formation is now 
on an exciting, rapidly developing stage. 
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